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UNIFORM UPPER BOUNDS FOR THE CYCLICITY OF THE 
ZERO SOLUTION OF THE ABEL DIFFERENTIAL EQUATION 


DMITRY BATENKOV AND GAL BINYAMINI 


Abstract. Given two polynomials P, q we consider the following question: 
“how large can the index of the first non-zero moment m*. = f b P k q be, 
assuming the sequence is not identically zero?”. The answer K to this question 
is known as the moment Bautin index, and we provide the first general upper 
bound: K ^ 2 + degg + 3(degP — l) 2 . The proof is based on qualitative 
analysis of linear ODEs, applied to Cauchy-type integrals of certain algebraic 
functions. 

The moment Bautin index plays an important role in the study of bifurca¬ 
tions of periodic solution in the polynomial Abel equation y' = py 2 + eqy 3 for 
p, q polynomials and 1. In particular, our result implies that for p satisfy¬ 
ing a well-known generic condition, the number of periodic solutions near the 
zero solution does not exceed 5 + deg q + 3 deg 2 p. This is the first such bound 
depending solely on the degrees of the Abel equation. 


1. Introduction 

1.1. Polynomial moments and their Bautin index. Throughout this paper 
P, Q G C[z] will denote a pair of polynomials, and p, q their respective derivatives. 
We will denote the degrees of P, Q (resp. p, q) by dp, (Iq (resp. d q ). We also fix 
two points a, b G C. 

Two related types of moment sequences corresponding to this data have been con¬ 
sidered in the literature, 

r b 

m k = m k (P,Q):= P k (z)Q(z)p(z)dz, 

J a 
fb 

m k = mk(P, q) '■= / P k {z)q(z)dz, 

J a 

These two sequences are related through a simple formula due to [Tjl (see 1J5] for 
details), and may be viewed essentially as different normalizations. The former nor¬ 
malization is particularly convenient for the study of moment generating functions, 
as explained in (J21 The latter appears more directly in the study of perturbations 
of the Abel equation, as explained in m 

Definition 1 . We define the vanishing index N(P, Q,a,b) to be the first index k 
such that m^(P, Q ) ^0, or oo if no such k exists. We define the moment Bautin 


fc = 0,1,2,... 

( 1 ) 

ft = 0,1,2,... 

( 2 ) 
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index N(dp,d,Q,a,b) to be 

N(d P , d,Q, a, b) := sup N(P,Q,a,b) + 1, (3) 

deg Q^cLq ,deg P^.dp 
N (P,Q,a,b)< oo 

i.e. the least k € N with the property that N(P , Q, a,b) ^ k implies N(P , Q , a, 6) = 
oo for any P, Q with deg P ^ dp and deg Q ^ g?q . 

We define N(P,q,a,b) and N(dp,d q ,a,b) analogously. 

Remark 2. Tfte moments mk(P,Q) are polynomials in the coefficients of P,Q. 
Let R denote the ring of polynomials in these coefficients and R C R denote the 
ideal by mo,..., m*,. Then N(dp , dg, a, 6) defined above is the first index for which 
the chain {\/Ifc}fceN stabilizes. In particular, from noetherianity it follows that this 
index is well-defined (finite). An analogous remark holds for N(dp,d q ,a,b). 

The moment Bautin index has been studied in various special cases, motivated 
primarily by its relation to perturbations of the Abel equation (see 1)1.21 for an 
overview). Bounds have been obtained in various special cases, including the cases 
dp = 2,3. We refer the reader to 011] and references therein for details. However, 
to our knowledge no general bound has been available. Our main result is the 
following general bound for the moment Bautin index. 

Theorem 1. For any dp,dQ £ N we have 

N(dp, d Q ,a, b) ^ d Q + 3(dp — l) 2 . (4) 

Similarly, for any dp,d q we have 

N(d P ,d q ,a,b) ^2 + d q +3(d P -l) 2 . (5) 

It is shown in 1[2] (following jTO]) that the second bound in Theorem [l] follows 
immediately from the first. Our approach to the proof of the first bound is based 
on the following two observations: 

(1) The vanishing index N(P,Q,a,b) is essentially equivalent to the order of 
the zero at t = oo of the moment generating function H (t) for the moment 
sequence {m*,}. It turns out [10] that H{t) admits an analytic expression 
as a Cauchy type integral for the algebraic function Q{P~ 1 {z)). 

(2) The Cauchy type integral above satisfies a (non-homogeneous) linear dif¬ 
ferential equation of Fuchsian type [7] . 

The problem of estimating N(P,Q,a,b) is thus reduced to the study of the order 
of zero at t = oo of solutions of certain Fuchsian differential equations. A detailed 
analysis of the Fuchsian differential operator involved, and elementary considera¬ 
tions concerning its monodromy, allow us to give an a-priori upper bound for this 
order of zero, thus proving Theorem [T] 


1.2. Perturbations of the Abel equation. The classical Hilbert’s 16th problem 
asks for bounding the number of limit cycles, i.e. isolated closed trajectories, of the 
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polynomial vector field 


— = -y + F{x,y), 
^ =x + G{x,y). 


( 6 ) 


The closely related Poincare Center-Focus Problem asks for explicit conditions on 
the polynomials F, G in order for the system (JH]) to have a center. These problems 
remain widely open, although during the years many partial results have been 
obtained (see j 6 ] for an exposition). 

An alternative context for the study of the problems above is provided by the Abel 
differential equation, 


y' = p( x )y 2 + q{x)y 3 , x e [a, b] c R, 


( 7 ) 


where p, q can be polynomials, trigonometric polynomials or even analytic functions 
[ 8 [. A periodic solution in this context corresponds to solution y(x) satisfying 
y(a) = y(b), and a center corresponds to an Abel equation where every solution 
with a sufficiently small initial condition is periodic. The Abel equation analogue 
of the Hilbert 16th problem, known as the Smale-Pugh problem, is to bound the 
number of periodic solutions of © in terms of the degrees of p and q. It is generally 
believed that some (but not all) of the essential difficulties in the study of © can 
be observed in 0 , even when one restricts to the case of polynomial coefficients. 
On the other hand, the polynomial Abel equation allows for several important 
technical simplifications, and significant progress has been achieved for the Center- 
Focus in this context using tools from polynomial composition algebra and algebraic 
geometry mm- 

The Smale-Pugh problem for the polynomial Abel equation remains open. Its 
infinitesimal version, first suggested in [ 5 ], is as follows: 


Problem 1. How many periodic solutions can a small perturbation 

y' = P{x)y 2 + eq(x)y 3 , x € [a,b\ (8) 

of the “integrable” equation y' = p{x)y 2 have? 


This is an Abel equation analog of the “Infinitesimal Hilbert 16th problem” for 
which an explicit bound was obtained in [Tj. Following [5], in this paper we focus 
our attention on the periodic solutions bifurcating from the zero solution of ©. 

The unperturbed equation (e = 0) is a center if and only if p(x)dx = 0. Thus we 
may choose the primitive P such that P{a) = P(b) = 0. As in the classical case, 
the study of the bifurcation of periodic solutions as well as the center conditions 
for the perturbation © begins with the study of the first variation of the Poincar 
map. 

For technical reasons it is customary to consider the “reverse” map from time x = b 
to time x = a. Namely, let G{y) : (C, 0) —> (C, 0) denote the germ of the analytic 
map assigning to each initial condition yb the value G(yb) = 77 (a), where 77 is a 
solution of © satisfying 77 ( 6 ) = 7/7,. We may view G as a germ of an analytic 
function in the coefficients of the polynomials p , q and e as well. Fixed points of G 
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correspond to periodic solutions, and the identical vanishing of G{y) corresponds 
to a center. An explicit computation [5] Proposition 4.1] gives the expansion 

q{x) 




1 - yP{x) 


d X = Y1 


m k y 


k+3 


(9) 


k -0 


As in the classical study of perturbation of Hamiltonian planar systems, it follows 
from this variational computation that the number of periodic solutions bifurcating 
from the zero solution of (HJ) is bounded by the order of zero of the right hand side, 
i.e. 1V(P, q, a, b) + 3, assuming that this number is finite. On the other hand, if the 
first variation vanishes identically then one must in general consider higher order 
variations in e, further complicating the study of bifurcating periodic solutions. 

A surprising feature of the Abel equation ([ 8 ]) is that for many polynomials p, the 
vanishing of the first variation @ automatically implies the identical vanishing of 
the Poincar map. Toward this end we recall the following definition. 

Definition 3 ([3]). The polynomials P , Q are said to satisfy the composition con¬ 
dition (PCC) on [a, b] if there exists a polynomial W(x) with W(a ) = W(b), and 
polynomials P, Q such that 

P( x) = P(W(x)), Q(x) = Q{W{x)). 


A polynomial P is called “definite” (w.r.t a,b), if for any polynomial Q, vanishing 
of all the moments fh k (P, q) implies PCC for P, Q. 


Definite polynomials are ubiquitous. In the deep works mm all counter-examples 
have been classified and shown to admit a rigid algebraic structure. 

Whenever the polynomials P,Q satisfy the PCC, the corresponding Abel equa¬ 
tion © automatically admits a center, as can be seen by a simple change of variable 
argument. We thus see that for a definite polynomial P, the vanishing of all mo¬ 
ments mfe(P, q) implies the identical vanishing of the Poincar map G(y). Therefore, 
in a sense the bifurcation of periodic solutions in ((SJ) is fully controlled by the first 
variation ©• More formally, the following holds. 

Theorem 2 ([3]). Let P be a definite polynomial, and fix the parameters a, b, d q . 
Then for any ||g|| <C 1 with deg q ^ d q , the number of periodic solutions of 0 with 
|p(a)| -C 1 is at most N(dp , d q , a, b) + 3. 

As a Corollary of Theorem [T| we therefore have the following first general estimate 
for the number of limit cycles near the zero solution for an Abel equation 0 with 
||<?|| small. 

Corollary 4. Under the conditions of Theorem [H the number of periodic solutions 
is bounded by 5 + d q + 3 dp. 

1.3. Organization of the paper. In >J2] we introduce moment generating func¬ 
tions for the two moment sequences {m k }, {fh k } which turn out to be Cauchy-type 
integrals. We thus reduced the study of the corresponding vanishing indices to the 
study of the order of zero of these generating functions at infinity. In ]J5] we give a 
slightly generalized version of the result of [7] which states that if a function g(z) 
satisfies a linear ODE Lg = 0 then the corresponding Cauchy-type integral I(t ) 



BOUNDS FOR CYCLICITY OF ZERO SOLUTION OF ABEL EQUATION 5 

satisfies a non-homogeneous linear ODE LI = R, where R is a rational function 
of known degree. Subsequently, in SJH we explicitly derive the corresponding non- 
homogeneous ODE for the moment generating functions. Finally in )}5]we produce 
estimates for the order of zero the moment generating function at infinity using 
qualitative methods of linear ODEs. 


2. Polynomial moments and generating functions 


Recall the notations of EH We introduce moment generating functions with the 
corresponding integral expression for the sequences {mk},{fhk} as follows: 


00 

H(t) = Y, m k t- {k+1) 

k =0 


( 10 ) 

00 

H{t) = Y fbikt~ {k+1) 

k=0 


( 11 ) 

Clearly, 



ordoo H (t) = N(P, Q, a, b) + 1 

ordoo H{t) = N{P,q,a,b) + 1. 

( 12 ) 


In particular, we have the following. 


Proposition 5. We have 

N(d P , d Q ,a, b) = sup ord t=oa H(t). (13) 

*nt)m 

where the supremum is taken over all pairs P, Q with respective degrees bounded by 
dp,dQ and H(t) denotes the corresponding moment generating function. 

It turns out that H(t) and H(t) are related by a simple formula, which implies 
in particular that the study of their orders of vanishing at t = oo are essentially 
the same [TO] Claim, p.40]. We repeat the argument of [ 10 ] in order to obtain an 
explicit description of relation between these orders. 

Lemma 6. The condition fho = 0 is equivalent to Q{a ) = Q{b). Moreover, under 
this condition we have fhk+i = — (fc + l)nik for k £ N. In particular, we have 

N(P, q, a, b) ^ N(P, Q, a, b) + 1. (14) 


Proof. Derivating under the integral sign we have 


d Hjt) 
d t 



Q(z)p(z) 

( t-P{z)Y 
Q{z ) -,b 


t-P(z) ia 


d * = -/>(d^) 

[ b q(z) , = Q{a) 
J a t-P(z) t — P(a) 


Q(b ) 

t-P(b) 


+ H{t) 


(15) 


Comparing the t 1 coefficient we see that mo = 0 if and only if Q(a) = Q(b), and 
under this condition rhk+i = — (k + 1 )mk as claimed. □ 


The moment generating function (fllll) has the form of a Cauchy integral. Indeed, 
choose the curve of integration 7 ' from a to 6 in (flUl) to be some smooth curve 
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avoiding the critical values of P(z) (except perhaps at the endpoints). Then setting 
7 = P( 7') and substituting w = P{z) in (TTT)]) we obtain 


H(t) = 


Q(P^H) dw 

t — w 


where P 1 (w) denotes the branch of P 1 lifting 7 to 7'. 


(16) 


3. Cauchy-type integrals and linear differential operators 

Let £ be a scalar differential operator, 

L = c r {z)d r H-+co(z), co,...,c r eC[4 (17) 

Let 7 C C be a smooth curve, and assume that 7 does not pass through the singular 
points of L, except perhaps at its endpoints. Finally let g be a solution of Lg = 0 
defined on 7 , and assume further that g is bounded on 7 (including at the possibly 
singular endpoints). We denote by p+,p- the endpoints of 7 . 

Then we define the Cauchy-type integral 

m = f^ m 

It is classically known that I[t) is a holomorphic functions defined on € \ 7 , and 
moreover that the boundary values I + and I~ of I(t) on 7 from above and below 
respectively satisfy I + — I~ = g | 7 . Moreover /(f) can be analytically continued 
along any path avoiding the endpoints of 7 . 

Kisunko [7j proved the following (under the extra mild assumption that g is holo¬ 
morphic at the endpoints of 7 ). 

Proposition 7. We have Lift) = Rft) where R(t) is a rational function having 
poles of order at most r at p+,p- and no other poles on C. 

Sketch of proof. By the classical properties of I ft) mentioned above, Lift) is a 
(possibly multivalued) analytic function on C \ {p+,p_} with ramifications p_(_,p_, 
and the difference between the two branches near the branch cut at 7 is g. But since 
Lg = 0, the boundary values of LI + and LI~ agree, so LI is in fact a univalued 
holomorphic function defined on C\ {p+,p_}. We will show that it has poles of 
order at most r at p+,p- and at most a pole at 00 . 

Since <?| 7 is bounded, we may derive under the integral sign and write 

<“> 

We now show that Lift ) admits polynomial growth of order at most r at p+ (and 
the same arguments work for p~). It is enough to consider each of the integrals 
in m separately. Moreover, we may assume that 7 is a small piece of a smooth 
curve near p + (because the integral over the rest of 7 is analytic at p + ). Choose a 
coordinate system where p+ = 0. 

Let M denote an upper bound for \g(z)\ on 7 . Let t be a point in a punctured 
neighborhood of p+. Since g(z) admits analytic continuation along any curve in 
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the punctured neighborhood, may deform 7 without changing LI(t) so that for 
some positive constants C , D independent of t, 

( 1 ) For every z £ 7 , we have \z — t\ ^ C\t\ and also \z — t\ > C\z\. 

(2) On 7 we have |d^| < Dd\z\. 

(3) Write g = 71+72 where 71 is the part of 7 which lies in {z : \z\ < 2 1} and 
72 is the rest. Then the length of 71 is at most D\t\, and the length of 72 
is at most D. 


We now estimate 




< length( 7 l )M(C|t |)- fe - 1 
proving the claim. 


72 {c\z\Y 
MDC k+1 , 


I — k 


< 0 ( 1 * 1 "*) ( 20 ) 


J2|t| 


Finally, it is easy to see that I it) and its derivatives have a zero at t = 00 , and 
since the coefficients of L are polynomial it follows that I(t) has at most a pole at 
00 as well. □ 


4. A DIFFERENTIAL OPERATOR FOR Q(P X ) 

Let V denote the linear space spanned by the dp branches of the algebraic function 
g{z) := Q{P~ 1 {z)). We denote r := dimt^, and note that r may be strictly smaller 
than dp. Denote by pi,... ,p s the critical values of P. 


4.1. The operator £. By a theorem of Riemann [BJ Theorem 19.7], there exists 
a linear r-th order differential operator £, with polynomial coefficients, 

£ = c r (z)d r H- \~c 0 (z), c 0 ,..., c r £ C[z] (21) 

whose solution space coincides with V. Moreover, L is uniquely determined by the 
requirement that c r ,..., Co do not share a non-trivial common factor. We recall the 
construction of L. 


Recall that the Wronskian W(/i,..., f n ) of a tuple of functions is defined to be 

( h ••• In \ 


:= det 


dfi 

V^-Vi 


df, 

d n ~ l fn) 


( 22 ) 


Now let gi,... ,g r denote r branches of g{z) which span V . Then clearly for any 
/ £ V we have W(gi ,..., g r , f) =0. We define the operator L given by 

r —1 


£(/) = 


W{gi,...,g r ,f) 

W r 


= [d r + ^2 c k {z)d k ] / where c k = 


k—0 


Wfc(ffi,...,g r ) 
W 7 r {gi , • • ■ > gr) 


( 23 ) 

where Wi are the minors obtained when expanding the Wronskian W (g 1 ,..., g r , f) 
along the last column. If the monodromy of g along a closed curve 7 induces the 
linear automorphism Ad . T : V —> V then the corresponding monodromy along 7 of 
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each Wk is given by multiplication by det M 7 . In particular, the coefficients Ck are 
uni valued functions. 

4.2. The divisors [Wk]- Let k = 0 , ...,r and zo £ CP. Choose any local rep¬ 
resentative of the functions g i,... ,g r . Since these functions have at most a finite 
ramification and moderate growth at zq, we may expand 

OO 

w k = J2 akj(z-z 0 ) ilg (24) 

j— — N 

where q and N are some natural numbers. Suppose that (ik,j a is the first non-zero 
coefficient among the akj ■ Then we say that the fractional order of Wk at Zq is 
ord , 0 W k '■= jo/q. This notion is well-defined: indeed, the monodromy of Wk along 
any curve is given by multiplication by a non-zero constant and hence does not 
change the order. We define the fractional divisor \Wk\ of W to be 

[W k ] := ^2 ord z W k {z)[z]. (25) 

zeCP 

This sum is locally-finite, and hence finite. Moreover it is clear that [cife] = \Wk] — 
[W r \. In particular, Ck admits finitely many singularities of finite order. Since we 
have already seen that Ck is univalued, it is in fact a rational function. 

We can also write the divisor [IT fc ] in terms of residues. Indeed, since the mon¬ 
odromy of Wk along any curve is given by multiplication by a constant, the one- 
form d LogWfc is a univalued one-form. It is easy to verify in local coordinates that 
it in fact has only finitely many poles, all of first order, and 


[W k ] = V Res(dLogtI4)[2]. (26) 

z ' 2 
zeCP 

For any divisor D = n* [zf\ we denote D Zi = m and 

degD = D + = ^ n t [zi\, D~ = - ^ n z [zi]. (27) 

In particular, it follows from the above that deg[114] = 0. 

4.3. An estimate for deg[I14] + . Our next goal is to estimate deg[I14] + . Since 
[IT r ] is principal, it will suffice to estimate deg[IT r ] _ . Recall that W r = W(g \,..., g r ) 
where gk = Ql y P^ 1 {z)) and Pff 1 {z),...,Pp 1 {z) denote r different branches of 
P~ 1 (z). If 2 £ C is not a critical value of P then these functions are all holo- 
morphic around 2 , and hence [W r ] z is non-negative. 

Let the critical value p t have exactly mt < dp preimages, and write bi := dp — m.i 
for the number of critical points (counted with multiplicities) over pi. Then at most 
2 bi of the branches gk may be ramified at Pi. We expand the determinant defining 
W r and note that: 

• since gk is bounded, its order is non-negative; 

• differentiation can decrease the order by at most 1; 

• differentiation cannot decrease the order below zero for holomorphic gk■ 
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We thus conclude that 


ord Pi W r > (—r + 1 ) H-+ (—r + u), where v = min(r, 2 bi). 

Since b\ + • • • + b s = dp — 1, it is not hard to see that the maximal value for the 
following sum is obtained when bi = 1 for i = 1 ,..., s, and in any case 

S 

Y ord Pi W r > - (2 d P - 3) (d P - 1). (28) 

i=l 


It remains to estimate the order of W r at oo. Choose a coordinate w around oo such 
that P(w) = w~ dp . Then any branch of Q(P~ 1 (w)) has the Puiseux expansion 

Q(p- 1 (w)) = Q(w- 1/dp ) = aw~ d Q /dp +■■■ , «/0 (29) 


where • ■ ■ denote higher order terms. Moreover, the derivative d z = —w 2 d w in¬ 
creases the order of zero at w = 0 by at least one. Expanding the determinant 
defining W r we see that 


_ r r t r do fir — 1 ) 

or doc W r ^ "t- 2 -' 


(30) 


In conclusion, we have 


degfWr]' 1 " = deg[W r ] < + (2dp - 3)(dp - 1) - . (31) 

ap A 


4.4. An estimate for degc r . We wish to derive an estimate for the number of 
singularities of L, or more specifically for degc r . By definition, c r is a polynomial 
and [c r ]+ is the least common upper bound for [co][c r _i] _ in (l23l) . Recall 
that [c fe ] = [W k \ - [W r ], 

We hrst note that £ is a Fuchsian operator. Indeed, since the solutions of £>, being 
algebraic functions, have moderate growth at each singularity, this follows from a 
theorem of Fuchs [BJ Theorem 19.20]. Thus by definition the order of [c r ] at any 
point p £ C cannot exceed r. We will apply this to the points pi,... ,p s . 

Let now z £ C and 2 ^ {pi,... ,p s }. Then the branches g±,... ,g r are holomorphic 
at z, so z is not a point of [Wo] - ,.. ., [W r - 1 ] - . In other words, z can only be a 
point of [co] - ,..., [c r _i] - if it comes from [W r ] + . Thus we see that 

S 

[ Cr ]+ < y r ip>} + [^] + ( 32 ) 

i=1 

Using (1311) and noting that s ^ dp — 1 and r ^ dp, we have the following Proposi¬ 
tion. 


Proposition 8. The following estimate holds, 


deg c r = deg[c r ] + < + 3 (d P - l ) 2 - ———. 

Up 


2 


( 33 ) 
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5. Estimate for the order of H{t) at infinity 

Recall from (£3 that H(t) denotes the moment generating function (1101) , which can 
be represented (around t = oo) as a Cauchy-type integral (fl6l) of the algebraic 
function g(z) = Q(P~ 1 (z)). Recall from §4] that V denotes the linear span of the 
branches of g(z) with r := dimU and L the differential operator (1211) satisfying 
V = ker L. 

Proposition 9. If H(t) ^ 0 then LH{t) ^ 0. 

Proof. Assume that LH(t) = 0. Then H(t) G V. Moreover, H(t) is holomorpliic 
at t = oo, and in particular it is invariant under the monodromy around infinity 
Moo and hence also under the operator 

^ dp — 1 

Too-.V^V, ^ := j- Y, M ( 34 ) 

P fe=o 

Recall that g{z) = Q(P~ 1 (z)) and P _1 (z) has cyclic monodromy at oo. It follows 
that the image of Too is one-dimensional and spanned by 

lmToo=C{S}, S(t):= ^ Q(w). (35) 

w:P(w)=t 

Moreover, S(t) is a polynomial: for instance, it is has no poles on C and moderate 
growth at oo. We conclude that H(t ) is a polynomial. Finally, H(t) has a zero 
at t = oo by definition, and since it is also a polynomial it follows that H(t ) = 0, 
contradicting the hypothesis. □ 


Let D = td t denote the Euler operator, and recall that it also gives the Euler 
operator at t = oo (up to a sign). The following proposition describes the behavior 
of L around infinity. 


Proposition 10. We may write 

L(t) = u(t)L 2 :=D r + c r -iD r ~ 1 + --- + c 0 , (36) 


where c r _i,...,co gwt rational functions, holomorphic at t = 
rational function satisfying 


ordoo u ^ 


dQr 

dp 


+ 3(d P — l) 2 


r(r + 1) 
2 


oo, and u(t) is a 


(37) 


Proof. The existence of an expression (1361) is a direct consequence of the fact that 
L is a Fuchsian operator at t = oo (see [ 6 ;, Proposition 19.18]). Using Proposition [ 8 ] 


we have 


ordoo u = r — deg c r ^ 


~T~ + 3(dp — l ) 2 

dp 


r(r + 1) 

2 


(38) 


as claimed. 


□ 


Finally we have the following estimate. 
Lemma 11. If H(t) ^ 0 then 
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Proof. Using Proposition 0] we have LH(t) = R(t), where R(t) has at most two 
poles of order r in C. Moreover, by Proposition [9] R(t) is non-zero. It follows that 
ordoo R(t) < 2 r. Using Proposition [TO] we have 


ordoo(££f(f)) = ord^ R(t) - ordoo u{t) 

<2 r+^: + 3(<lp-l f-±±A 


It remains only to note that the application of L cannot decrease the order of zero, 
and the claim follows. □ 


Finally we complete the proof of our main result. 

Proof of Theorem 0 If H{t) ^ 0 then by Lemma ITT! 

ordoo H(t) < dQ + 3(d P - l) 2 , (40) 

and the claim for N(dp, d.Q,a, b) follows by Proposition[5] The claim for N(dp, d q , a , b) 
then follows from Lemma [Gl noting that cIq = d q + 1. □ 
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